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ABSTRACT 

We develop population models of high-mass X-ray binaries (HMXBs) formed after bursts of star 
formation and we investigate the effect of electron-capture supernovae (ECS) of massive ONeMg white 
dwarfs and the hypothesis that ECS events are associated with typically low supernova kicks imparted 
to the nascent neutron stars. We identify an interesting ECS bump in the time evolution of HMXB 
numbers; this bump is caused by significantly increased production of wind-fed HMXBs 20-60 Myr post 
starburst. The amplitude and age extent of the ECS bump depend on the strength of ECS kicks and the 
mass range of ECS progenitors. We also find that ECS-HMXBs form through a specific evolutionary 
channel that is expected to lead to binaries with Be donors in wide orbits. These characteristics, 
along with their sensitivity to ECS properties, provide us with an intriguing opportunity to probe 
ECS physics and progenitors through studies of starbursts of different ages. Specifically, the case 
of the Small Magellanic Cloud, with a significant observed population of Be HMXBs and starburst 
activity 30-60 Myr ago, arises as a promising laboratory for understanding the role of ECS in neutron 
star formation. 

Subject headings: galaxies: starburst — stars: emission-line, Be — stars: neutron — supernovae: 
general — X-rays: binaries 



1. INTRODUCTION AND BACKGROUND 

Recently, electron-capture supernovae (ECS) of mas- 
sive ONeMg white dwarfs have attracted renewed atten- 
tion as an important mechanism of neutron star (NS) 
formation. ECS events probably operate at the low end 
of the NS progenitor mass range, when a strongly de- 
generate ONeMg core forms. If the progenitor exceeds 
a critical core mass of 1.38 M , high pressure and den- 
sity lead to electron captures onto 24 Mg and 20 Ne. The 
resulting decrease in electron degeneracy pressure causes 
a super nova (SN) co l lapse before an iron core is ever 
formed (jMivaii et all [l980t iNomotol [l98l H98l . Nu- 
merical simulations indicate that ECS events are under- 
energetic compared to iron-core supernovae (ICC-SN 1 ) 
by at least an order of magnitude, However they are, 
at the same time, relatively fast, neutrino-driven, de- 
layed explosions. The short time available for the prop- 
agation of the shock is expected to prevent a strong 
build-up of core asymmetries, leading to typical kicks 
of < 100 km s" 1 , lower than those in ICC-SN events 
(jDessart et alj2006tlKitaura et al.ll2006l Janka 2007, pri- 
vate communication). 

The mass range of ECS progen itors is highly uncer- 
tain. While INomotol (|1984 11987ft predicted an initial 
mass range of 8-10 M Q at solar metallicity, recent simu- 
lations suggest mu ch stricter ECS boundar ies for single 
stars. For example. iPoelarends et all (|2008f ) find a single 
star mass range of 9-9.25 M©, which accounts fo r only 
~4% of SNe for a given population. iSiessI (|2007l ) finds 

1 Here we use the term ICC to refer to both iron-core collapse 
supernovae, as well as direct collapse black holes (BHs) which form 
without a SN explosion (i.e, all non-ECS compact object formation 
events). We use ICC-SN to refer to only those iron-core collapse 
systems which result in an energetic supernova explosion 



a metallicity-dependent mass range which stays consis- 
tently narrow (a width of <1.5Mq), and could vanish 
altogether depending on the metallicity dependence of 
t he stellar winds of ev o lved s tars . 

iPodsiadlowski et all (2004) proposed that binary evo- 
lution can result in a much wider ECS initial mass 
range of 8-11 M Q in the following way. Stripping of 
the hydrogen-rich stellar envelope due to mass trans- 
fer and/or common envelope evolution aborts the sec- 
ond dredge-up phase, normally seen in single stars with 
their hydrogen envelope intact. Consequently the mass 
of the helium core is preserved, leading to relatively mas- 
sive ONeMg cores that experience electron-capture col- 
lapse. They find that, in the closest binaries, even the 
first dredge up phase may be aborted, enabling ECS ac- 
tivity in the even lower mass range of 6-8 M Q . 

Given the uncertainties involved, probing the ECS 
physical mechanism and its relevance to NS formation 
compared to ICC-SN events is a challenge. If ECS events 
are indeed associated wi th lower kicks (|Dessart et al.1 
120061 : iKitaura et al.ll2006l ). then close binaries contain- 
ing a NS are good candidate probes as their formation 
rates and prope rties signific a ntly d epend on SN kicks 
and mass loss. iPfahl et al.l (|2002T ) highlighted a sub- 
class of high-mass X-ray binaries (HMXBs) with low 
eccentricities that c a n be expla ined by low NS kicks . 
IPodsiadlowski et all (|2004h and Ivan den Heuvell (|2004l 
l2007h postulated that the ECS mechanism could be re- 
sponsible for these low kicks, although a very wide ECS 
progenitor mass range may be required to account for 
the observed HMXBs. 

In this paper, we present results that uncover an alter- 
native probe for ECS events through studies of HMXB 
formation in starbursts of varying ages. We find that 
ECS NS formation associated with low kicks (lower 
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than ~50kms _1 ) dramatically increases the formation 
efficiency of Be X-ray binaries in starburst ages of ~ 
20 — 60Myr. In what follows, we discuss this effect in 
detail and analyze the reasons for this unique depen- 
dence, with the hope that observations of starbursts can 
be used to assess ECS kick velocities and progenitor mass 
ranges. We note that, in the Small Magellenic Cloud 
(SMC), dominated by ~ 40Myr-old stellar populations, 
a syste matic overabundance of Be XRBs has b een ob- 
served (jHaberl fc Sasakill2000t iMaiid et al.ll2004h poten- 
tially making the SMC a unique laboratory for ECS NS 
formation. 

2. SIMULATION CODE AND MODELS 

For our modeling of HMXB formation and evolution 
in starbursts we employ a sophisticated population syn- 
thesis code, Star - Track , described in extensive detail in 
iBelczvnski et al.l ((2008), which assumes that stellar dy- 
namical interactions are not significant compared to the 
effects of binary evolution for the star formation con- 
ditions under consideration. We note that the param- 
eter space of Monte Carlo population synthesis is very 
large, and thus a full exploration is not possible. In- 
st ead, we consider a de fault model (described in detail 
in IBelczvnski et all (|2008D ) and we vary some binary evo- 
lution and ECS-related parameters with the purpose of 
highlighting their effects on the proposed ECS probe. 
Here we briefly summarize the main assumptions and 
parameters relevant to this study. 

We employ a delta function star formation episode at 
solar metallicity adopting: (1) a Salpeter (M -2 35 ) initial 
mass function with primary masses above 4 Mq and sec- 
ondary masses above 0.15 M Q ; (2) a flat mass ratio dis- 
tribution; (3) a distribution of initial binary separations 
that is flat in the logarithm with an upper limit of 10 5 Rq 
and a lower limit such that th e primary s tar initially fills 
at most half of its Roche Lobe ((Abt 1983), and (4) a ther - 
mal distribution for initial eccentricities (Hcggie 1975). 
We set the maximum NS mass at 2.5 M© and draw 
natal kicks for "standard" ICC-SN events from a sin- 
gle Maxwellian kick distribution with mean 265 km s _1 



(Hobbs et alj [2005) . Kicks ar e potentially associ ated 
with BH formation as well (see iFragos et al.l t|2009h for 
the strongest evidence at present for a BH kick). For 
BHs formed through SNe explosions and subsequent fall- 
back of material we multiply the normal Maxwellian kick 
by the fraction of the SN ejecta which is ultimately lost 
from the system. Usually BH formation at higher masses 
through direct collapse is assumed to be perfectly sym- 
metric; here we adopt a small kick (10% of the NS kick), 
for reasons discussed in detail in Section Q] 

To obtain statistically significant results, we sample at 
least 10 6 initial binaries for each set of input parame- 
ters and we evolve the systems for 100 Myr. However, 
we stress that the models and XRB numbers presented 
here are not normalized to match any specific observed 
system, since we are purely interested in examining the 
relative XRB numbers as a function of starburst age. 

In order to compare our results with Chandra obser- 
vations, we must determine the X-ray luminosity, Lx, 
based on the calculated mass-transfer rate for both wind- 
fed and Roche-lobe-overflow systems. We apply an esti- 
mated correction for t he Chandra energy ban d as de- 
scribed and justified in IBelczvnski et ail (20081 ) (Section 
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Fig. 1.— Number of HMXBs with L x > 1 x 10 32 ergs _1 per 
10 6 initial binaries vs. starburst age for: (top) different Maxwellian 
sigma values for ECS natal kicks, and (bottom) different core mass 
ranges at the base of the AGB (BAGB) and no ECS kicks 

9.1). We analyze the HMXB population with Lx in 
excess of 10 32 erg s , appropriate for observations of 
nearby star-forming galaxies. 

An important problem in the population synthesis of 
binary stars is the treatment of common envelope (CE) 
phases. In our simulations, CE events are tr eated us- 
ing th e usual energy formalism described by IWebbinkl 
(1984) 2 , for stars which have established a clear core- 
envelope boundary. In this paper, we adopt a value 
for the CE efficiency of acE = 1- However, in our de- 
fault models, we assume that a CE phase involving a 
donor star on the main sequence ( MS) or hertzsprung gap 
(HG) leads to a binary merger (|Belczvnski et alj 120081 : 
iTaam fc Sandauistll2000D 

Following iHurlev eTatl HH), massive stars are as- 
sumed to explode in ECS events if the He core mass at the 
beginning of the asymptotic giant branch (AGB ) is be- 
tween 1.83-2.25 M (see also llvanova et al.l (j2007T )l. This 
choice, in effect, selects specific initial mass ranges which 
depend on the binary evolution of the progenitors. How- 
ever, given the uncertainties involved in selecting this 
mass range, we examine a number of He core mass ranges 
in what follows. We assume that ECS natal kicks follow a 
Maxwellian distribution with a smaller mean than ICC- 
SN kicks, thus we linearly scale down the "standard" 
ICC-SN Maxwellian distribution by varying factors. The 
mass of ECS-formed NS is assumed equal to 1.26 M w 
(lNomotolll9"87t iPodsiadlowski et al.112001 llvanova et all 
20071). although we not e that iDessart et alJ (12006ft . 



Kitaura et all (20061) and iPoelarends et all (|2008l ) find 



slightly higher end masses for single star ECS events 
(1.338 - 1.4 M ). 

3. THE ROLE OF ELECTRON-CAPTURE SUPERNOVAE IN 
HMXB PRODUCTION 

Examination of our simulation results in terms of the 
number of HMXBs as a function of starburst age reveals 

2 In this treatment of CE events, the energy needed for CE 
ejection is removed from the binary orbital energy, often resulting 
in evolved stellar systems with very tight orbits. 
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an intriguing role of ECS NS formation: ECS events 
can dominate HMXB formation between 20 and 60 Myr 
post-starburst, creating a clearly identified "bump" in 
the time evolution of HMXB numbers (see Figure [1]). At 
any given metallicity, the width and relative height of 
this ECS bump is primarily dependent on two parame- 
ters: the typical ECS natal kick magnitude and the mass 
range of ECS progenitors. 

We first examine the effect of ECS kicks (Figure [TJ 
top panel). We typically find two significant bursts of 
HMXB activity: (1) a sharp, narrow peak at 5— 10 Myr 
post-starburst dominated by wind accretion from mas- 
sive donors onto massive black holes formed in "stan- 
dard" ICC events. As a result, this peak is, of course, in- 
sensitive to assumptions about ECS kicks; (2) a broader 
ECS bump (at 20-60 Myr) dominated by wind accretion 
from lower mass donors onto NS. The sensitivity of this 
bump on ECS kicks is evident, as it effectively disappears 
and blends into the background HMXB population once 
typical ECS kicks significantly exceed ~ 50kms _1 . 

As one might naturally expect, the ECS bump ampli- 
tude and age spread also depends on the AGB He core 
mass range for ECS events (Figured! bottom panel). We 
find that more massive cores lead to earlier ECS bumps, 
as higher-mass cores select more massive stars that form 
NSs earlier. The effects of binary evolution change the 
mapping between He core masses and zero-age MS pro- 
genitor masses for ECS. The details of this mapping will 
depend on the specific stellar models adopted, but for 
the three mass ranges shown in Figure Q] (bottom panel) 
and the StarTrack models, the corresponding progeni- 
tor masses are found in the following ranges: 7.5-8.5 M Q 
and 9-11 M Q (solid), 11-14.5 M (dotted), and 9-9.5 M 
and 11. 5-1 2.5 Mq (dashed). The change in the number of 
bright HMXBs for different ECS mass regimes by nearly 
a factor of 10 conveys the potential of this bump for 
constraining the details of the ECS mechanism and its 
dependence on the physical properties of the collapsing 
star. 

Analysis of the evolutionary pathways of HMXBs re- 
veals that CE evolution strongly governs the parame- 
ter space of bright HMXB progenitors. CE events lead- 
ing to binary mergers greatly limit this parameter space 
by effectively imposing a minimum initial binary sepa- 
ration for avoiding CE mergers. Wider systems (with 
more evolved donors) may survive the CE event, acquir- 
ing tight post-CE orbits which contribute to HMXB pro- 
duction in two ways: (1) for CEs occurring prior to NS 
formation, tight orbits promote survival of the SN kick 
and (2) for CEs occurring after NS formation, tight orbits 
lead to stronger wind accretion and brighter HMXBs. 

As already mentioned, in our default sim u lation s 
we follow the finding s of iTaam fc Sandquisll ([2000) ; 
llvanova fc Taa m (2004) and assume that CEs driven by 
MS or HG donors always lead to binary mergers. To 
examine the robustness of the ECS bump, we consider 
models where CEs driven by HG donors are also treated 
with the energy formalism. In this case more HMXBs 
form, but the presence of the ECS bump remains robust 
(Figure [2]) . Careful examination of the results indicates 
that the HMXB numbers populating both peaks (narrow 
and ECS) remain roughly constant, and the details of 
the CE treatment affect only the number of background 
(outside the peaks) systems throughout post-starburst 
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Fig. 2. — Number of HMXBs with Lj >lx 10 d2 ergs -1 when 
CEs involving HG donors are always assumed to lead to binary 
mergers (solid curve) and whe n CEs with HG donors are treated 
with the usual Wcbbink (1984) energy formalism, which may lead 
to either mergers or tight post-CE binaries. 

evolution. 

We find that the resilience of both our BH-HMXB and 
ECS-HMXB peaks to CE events stems from their evo- 
lution through a particular channel and mass transfer 
history. Beginning with initial mass ratios close to unity, 
and relatively tight orbits, these systems go through mass 
transfer before BH/NS formation, but avoid a CE phase. 
Instead, they evolve through a mass ratio inversion which 
significantly expands their orbits through prolonged sta- 
ble mass transfer. This pathway is dominant due to three 
main effects related to the initial binary properties, CE 
outcomes, and kick distributions of ECS and massive BH. 

The logarithmically flat distribution of initial orbital 
separations tends to favor tight systems. These un- 
dergo Roche-lobe overflow when the primary is moder- 
ately evolved (late on the MS, on the HG, or early in 
the core He burning [CHeB] phase). Given the merg- 
ers associated with CE events, we find that the majority 
of the surviving binaries have comparable masses (with 
accretor mass at least half the donor mass), ensuring sta- 
bility of the mass transfer initiated by the original pri- 
mary and avoiding a CE merger. Stable mass transfer 
soon reverses the binary mass ratio and the orbits begin 
to expand. This expansion continues until the hydrogen 
envelope of the original primary is exhausted, creating 
systems which are both evolved and relatively loosely 
bound (orbital separations of ^500-2500 R©). Conse- 
quently, when the original primary explodes in a SN, bi- 
nary survival is strongly favored if SN natal kicks are ab- 
sent or low, which is the case for massive BH and for NS 
formed through ECS events. HMXBs in the ECS bump 
(Figure [l} have relatively wide orbits and reach signifi- 
cant X-ray luminosities only when the donors (original 
secondaries) evolve off the MS onto the HG and CHeB 
phases where stellar winds are sufficiently strong. Lastly 
we note that, because of the comparable masses of the 
initial binary components, the post-MS evolution of the 
donor follows soon after the SN explosion of the primary, 
thus dictating the timing of the two identified peaks. 
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Through this evolutionary sequence, ECS events (as- 
sumed to be associated with relatively low SN kicks) 
lead predominantly to the formation of wide HMXBs 
with evolved donors that have reached the HG or the 
CHcB phase. The progenitors of these ECS-HMXBs 
have avoided a CE phase before the formation of the com- 
pact object and instead go through stable mass trans- 
fer from the original primary. Some systems can fur- 
ther evolve through a CE phase and form tightly bound 
HMXBs with He-rich donors, the relative contribution of 
which depends strongly on the assumed CE efficiency. 

On the other hand, binaries that evolve through a 
CE phase before compact object formation end up with 
tight enough orbits to survive ICC-SN events with typ- 
ical natal kicks in excess of ~ 100 kms . These bi- 
naries tend to remain tight after the SN explosion and 
can become bright HMXBs even with unevolved, MS 
donors. These systems form a background of short period 
HMXBs through which the ECS bump can rise depend- 
ing on the how favorable ECS conditions are (progeni- 
tor mass range and natal kicks). Since the creation of 
post-CE HMXBs is not dependent on the evolution of 
the donor, most MS-HMXBs form immediately after the 
SN of the primary (< 25Myr) and their numbers expo- 
nentially decay thereafter due to subsequent Roche-lobe 
overflow by the MS donor and ultimately mergers. 

Lastly, we note that the formation of systems with 
wide orbital periods (greater than ~ 10 days) is highly 
unlikely to result from ICC-SN activity involving large 
SN natal kicks (i.e. NS from non-ECS events and low- 
mass BH from ICC supernovae and fallback). Instead 
this population is dominated by ECS NS and BH formed 
through direct collapse, given the small kicks adopted for 
such events. However, as shown for ECS-HMXB in Fig- 
ured] (top), even these small kicks are enough to disrupt 
the widest of HMXBs. 

4. DISCUSSION 

We have shown that, if NS formation through ECS 
events is associated with natal kicks smaller than 
~ 50kms _1 , production of relatively bright wind- fed 
HMXBs (Lx > 1 x 10 32 ergs _1 ) is significantly favored 
between 20 and 60Myr after delta function star forma- 
tion. We call this increased HMXB formation, which 
appears above a background HMXB population formed 
through standard ICC compact object formation, the 
ECS bump. The width in post-starburst age of this ECS 
bump and its amplitude relative to the rest of the ICC- 
HMXB population 3 depend primarily on two ECS fac- 
tors: the mass range of ECS progenitors and the typical 
magnitude of ECS natal kicks (Figure QJ. Moreover, we 
find that these ECS-HMXBs form through a specific evo- 
lutionary channel that avoids a CE phase before the SN 
and instead includes stable mass transfer from the pri- 
mary and a mass ratio inversion leading to orbital expan- 
sion. As a result, ECS-HMXBs have predominantly wide 
orbits (~500-2500R Q ) and evolved donors. These char- 
acteristics, along with their sensitivity to ECS properties, 
provide us with an intriguing opportunity to probe ECS 
physics and progenitors through studies of starbursts at 
different ages. 

3 including a narrow peak of HMXBs with massive BHs that 
appears earlier at ~ 5 Myr 




1 10 1 00 1 000 1 4 1 5 

Distance Traveled before detection (pc) 

Fig. 3. — Distance travelled by HMXBs with L x > 10 32 ergs -1 
detected between 30-60 Myr after star formation for ECS SN with 
natal kick velocities of 26.5 km s _1 . Distances are calculated ig- 
noring any deceleration effects from the underlying galactic gravi- 
tational potential. 



The SMC is a promising candidate for a study prob- 
ing ECS physics. Observations of the SMC have re- 
yealed a major burst of star formation 30-60 My r ago 
(|Maragoudaki et al.ll200lt lHarris fc Zaritskvi 120041 ). Re- 
cent observations have found a surprisingly large popula- 
tion of HMXB s with Be donors conc e ntrated in the SMC 
bar region (lHaberl fc Sasakil l2000t iMaiid et al] 12004 
lAntoniou et al.ll2009fh ~ 

There are four reasons to prefer the hypothesis that 
Be HMXBs are primarily formed through ECS creation 
channels. First, the population of HMXBs in the SMC 
has been found to have wide orbital separations and 
low eccentricities, which would require fine tuning of 
the initial orbital parameters if strong kicks were used 
(Pfah l~et al.ll2002T ). Second, we note the coincidence of 
age between the peak of ECS driven HMXB activity, and 
the observed age of star forming regions of the SMC bar. 
Third, the spatial distribution of SMC HMXBs is highly 
peaked in the bar region, with few HMXBs elsewhere in 
the galaxy. This spatial distribution is highly unlikely to 
result from high velocity natal kicks oc curing between 10 
and 4 Myr before the present. Lastly. iMcSwain fc Giesl 
(2005) have recently argued that Be stars form through 
spin up due to mass transfer in binaries that experi- 
ence mass ratio inversion and associated orbital expan- 
sion, which fits the evolution pathway found for the ECS 
HMXBs in our models. 

In order to estimate the overconcentration of ECS 
sources in the SMC bar, we select a value of 
cecs = 26.5 km s _1 , and examine the number of HMXB 
with Lx > 1 x 10 32 ergs -1 between 30 and 60 Myr after 
star formation, taking snapshots of the systems position 
every 0.1 Myr. In Figure[3l we plot the cumulative distri- 
bution function of three-dimensional distances traveled 
due to the SN natal kick of both ECS and ICC sources, 
normalized to their total number. We ignore any effects 
due to the galactic gravitational potential on the motion 
of HMXB sources. We note that while overall there are 
approximately twice as many ICC sources, only around 
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15% of ICC sources are found within 1 kpc of the loca- 
tion in which they were formed. Inside this region, ECS 
sources dominate by a factor of approximately 4. We 
note that the deep Chandra observations concentrate 
on the SMC bar region, while regions further outside are 
either not targeted by observations or their X-ray lumi- 
nosity sensitivity limit is much higher. Consequently a 
significant number of XRBs predicted by the models to 
have moved to distances beyond ~ 10 kpc could remain 
undetected. 

Our above findings strongly depend on our choice to set 
natal kicks for direct collapse BHs to be 10% as strong as 
our Maxwellian natal kick distribution. If wc set direct 
collapse natal kicks to be strictly zero, we find a large 
population of persistent sources in the SMC bar with 
luminosity Lx > 10 35 ergs _1 , containing CHeB donors 
and massive BH accretors, and having orbital periods on 
the order of thousands of days. As these systems are not 
only unobserved in the SMC, but throughout all HMXB 
populations, we enable small direct collapse kicks which 
split these loosely bound systems. 

Lastly, while the Be star phenomena is not understood 
from first principles, making it difficult to determine the 
absolute number of systems with Be-driven luminosities 
above Lx > 10 34 ergs _1 , we note that our simulations 
are consistent with the lack of observed supergiants with 
spherically driven winds. For each million initial binaries, 
we find approximately 7.5 such sources with luminosity 
Lx > lO^ergs^ 1 which have travelled less than 1 kpc. 
Using standard assumptions for initial mass function and 
star formation rate in the bar region, this is consistent 
with a Poisson distribution of supergiant systems. 

These results demonstrate the promise of using the 
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ECS process in order to correctly model the great over- 
abundance of Be-HMXB sources in the SMC. Because of 
the low natal kicks imparted to ECS systems, HMXBs 
formed through the ECS process are able to replicate 
the long period of Be-HMXB sources and high concen- 
tration of these systems in the SMC bar. While several 
models exist which allow for low natal kick velocities, 
we show that the ECS process should be favored due to 
its ability to create a dominant population of systems at 
times matching SMC observations. In addition, the ECS 
mechanism produces HMXBs through the same physi- 
cal p rocess suspected to pro duce a Be-HXMB popula- 
tion (jMcSwain fc Giesl feOOS). Detailed modeling of the 
HMXB population in the SMC, taking into account the 
duration, age and spatial distribution of star formation 
activity, the motion of the star forming clusters, and the 
lower metallicity of the SMC is beyond the scope of this 
paper, but is currently underway. We anticipate that 
such studies of the SMC and potentially other nearby 
galaxies with recent starburst activity will greatly con- 
tribute to efforts of understanding SNe and compact ob- 
ject formation. 
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